Effect of Symmetry Distortions on Photoelectron Selection Rules and Spectra on 

Bi 2 Sr 2 CaCu 2 8+5 
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We derive photoelectron selection rules along the glide plane in orthorhombic Bi2Sr2CaCu2 08+a 
(Bi2212). These selection rules explain the reversed intensity behavior of the shadow and the main 
band of the material as a natural consequence of the variating representation of the final state as a 
function of fey. Our one-step simulations strongly support the structural origin of the shadow band 
but we also introduce a scenario for detecting antiferromagnetic signatures in low doping. 

PACS numbers: 79.60.-i, 71.18.+y, 74.72. Hs 



Angle Resolved Ultraviolet Photoelectron Spec- 
troscopy (ARUPS) is probably the most important tool 
to study the electronic structure of complex materials 
like the high T c superconductors. It has been shown 
that matrix element effects must be taken into account 
in interpreting experimental data^£. Using the one-step 
model of ARUPS and optical matrix elements, light has 
been shed on the strong variation of intensities due to the 
matrix element effects £iai£ However, there are more strict 
rules dictating intensites in ARUPS, selection rules due 
to crystal symmetry. High resolution angular resolved 
photoemission measurements along high symmetry lines 
have been carried out recently to investigate the origin 
of the shadow bands in Bi2212£*£. Based on one-step 
computations, it was proposed that the shadow Fermi 
surface(FS) in Bi2212 has its origin in the orthorhombic 
(orth.) structure of the material. 6 In this study a need for 
selection rules emerged. In the orth. structure a mirror 
plane is broken into a glide plane. In this letter we give 
a group theoretical analysis for the selection rules along 
the glide plane and demonstrate the derived rules with 
one-step ARUPS simulations. We focus on Bi2212, but 
the derived selection rules may be applied to any cuprate 
superconductor with a similar structure. 
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FIG. 1: Sketch of the positions of atoms in Bi2212. (a) Tetr. 
Bi-0 layer, (b) orth. Bi-O layer, (c) Sketch of the Fermi 
surface. 



y axis is changed into a glide plane which is plotted with 
dash-dotted line in Fig. 



1(b) This glide plane corre- 



The space group of the commonly used body-centered 
tetr agonal (tetr.) structure for Bi2212 is number 139 
in the International Tables for Chrystallography. In x- 
ray diffraction experiments an orth. structure has been 
found. % In this structure atoms are moved from their tetr. 
positions and the symmetry group is changed. The non- 
symmorphic space group of the orth. structure is the 
group number 66. The geometric structure of the super- 
conducting cuprates may be described by a pile of planes 
of atoms. In the orth. structure the movements of the 
atoms from tetr. positions in the Cu-O-planes are small. 
Displacements are more clear in the Bi-O-plane, which is 
scetched in Fig, fl] Fig. 1(a) represents the tetr. Bi-O- 



sponds the glide-operation {<7j;|b/2}, a reflection (a x in 
the x — plane) followed by a translation by half the 
primitive lattice vector in the y axis. 

The model FS of Bi2212 (only one sheet is shown for 
clarity) is sketched in Fig. |l(c)| The figure shows the 
main band (FS of the tetr. structure) as thick circles 
and the shadow band as thin circles. In the figure the 
first tetr. Brillouin zone is approximately plotted with 
dashed line and the first orth. Brillouin zone with solid 
line. 2-dimensional fc|| -points of high symmetry T,X and 
Y are also marked. The direction from T to X that is 
parallel to the mirror plane in real space is denoted by 
E, and the direction from f to Y that is parallel to the 
glide plane by Am> These directions are commonly known 
as the nodal directions. 

The formula, based on the Fermi's golden rule, for the 
photocurrent by Feibelman and Eastman^ can be ma- 
nipulated into the form 



/(^A^cx/m^V^IA-PhA. 



(initial) 



plane and Fig. 11(b)! the orth. Bi-0 plane, where empty 



circles representing Bi-atoms and filled circles represent- 
ing O-atoms are displaced. In the tetr. structure there 
are mirror planes corresponding to dashed lines in Fig. 
|l(a)[ but in the orth. structure only the mirror plane 
parallel to x axis remains. The mirror plane parallel to 



{initial) 



A ■ p|^> 



(final) 
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where A is the vector potential of the incident photon 
and Gij is a hole Green function. Photoemission selec- 
tion rules arise from the symmetries of the initial (jpi) 
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and final (ipf) states and the dipole operator A • p in the 
matrix elements of Eq. ^ These selection rules are es- 
pecially relevant along the mirr or and glide planes, the 
directions £ and A in Fig. 1(c) In the following analy- 



sis the dipole approximation is used, but the derivations 
are exact beoynd the approximation if the wave vector 
of the radiation lies in the symmetry plane that is under 
discussion. 

In the mirror plane, since there is a well defined parity 
about the reflection in the plane, the selection rules can 
be derived by considering the parities of the components 
of the matrix element^ The transition is allowed if the 
entire dipole matrix element (ipf\A ■ p\ipi) is even. With 
respect to reflection the final state with its momentum 
lying in the mirror plane is even. The operator A • p is 
even(odd) if the polarization is parallel(perpendicular) to 
the plane, which can be explained by considering parities 
of the initial state ipi an( i its derivative A • pipi, where 
p = —ifiV. The initial state at the Fermi level is sim- 
ilar to a Bloch sum of copper d x 2_ y 2 orbitals 20 . With 
respect to reflection this orbital and the Bloch sum with 
its momentum lying in the plane is odd. Thus, for pho- 
tointensity to be nonzero the operator A ■ p must be odd 
and respectively polarization must be perpendicular to 
the mirror plane. This result is presented in combined 
form in table [D] The even(identity) representation is de- 
noted by £ and the odd representation by £'. 

Along the glide plane in the orth. structure the full 
machinery of the group theory has to be used. Pho- 
toemission selection rules due to glide plane have been 
considered previously by Pescia et alH A good account 
of a more general procedure has been given by Bassani 
et alM. The dipole matrix element |(A • p)( a ) \tpi v) ) 
will vanish if the irreducible representation(irrep) of the 
final state does not appear in the product of irreps of the 
initial state and the operator. The number of times it 
appears is given byl^ 

c(k, fi; a; k, v) — 

^ £ X (Al) ({i?k|f})*X (Q) (^)x (l/) ({^k|f}), (2) 



where hk is the order(number of elements {i?k|f}) of 
the relevant little group for a particular fe-point in the 
first Brillouin zone, x ({-Rk|f }) is the character of 
the little group irrep ({Rk\f }) for the initial state, 
({^k|f }) the character of the irrep of the final state 
and x^ (R) is the character of the small point group rep- 
resentation D^ a \Rk) of the dipole operator. The little 
group of a particular k point contains operations {i?|f} 
of the full space group that satisfy 



Rk = k + G, 



(3) 



i.e, leaving k unchanged. The small point group of k con- 
tains the rotational parts of the little group operations. 
In the glide plane determined by direction A in Fig. 



k z unchanged) of the k vector in the plane contains two 
operations, identity operation {-E|0}, and glide operation 
{(j^lb/2}. Both of the two operations form a class of their 
own and the eigenstates in the plane can be classified by 
two irreps. Character table for the little group is in table 

m 



TABLE I: Character table of the little group in plane A. S is 

-ikb/2 
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K|b/2} 
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S 


A' 


1 


-S 



To proceed it has to be found out how to assign sym- 
metry labels to particular final (and initial) states. The 
matching of the wave function at the surface of the sam- 
ple rules out some of the possible final states. The final 
state in ARUPS is a time-reversed LEED stateiS. The 
irrep of this state may be obtained from the plane wave 
state ij)k } = e l ( k/ ' r ) that arrives to the detector—. The 
state is transformed under the glide operation according 
to 



0{a x \b/2}ipk f (r) = tpk f (-x,y- 6/2, z) 
= e- l(fc « +G ^ /2 ?/> fc/ (r) = 6e- m ^ kf (r), 



(4) 



1(c) the little group (operations that keep fen = fc y k y and 



where k y lies in the first Brillouin zone and G is a recip- 
rocal lattice vector. The glide operation is represented by 
(5e~ m7r , where n is integer. The irrep of the final state is 
a function of the magnitude of the parallel component of 
its wavevector. It belongs to the irrep A in the l.,3.,5.,... 
repeated Brillouin zone, and to the representation A' in 
the 2., 4., 6.,... zone. This fundamental result for ARUPS 
and glide plane was first shown by Pescia et alH and has 
been exploited for Ni(100)-p(2x2)C±i and TiOiA 

Operator A • p belongs to representation A of the small 
point group if polarization vector lies in the glide plane 
and to representation A' if it is perpendicular to the plane 
(parallel to x-axis). Photoemission selection rules follow 
from Eq. [21 They are presented in table ^ In summary 
the initial state must belong to the same irrep as the final 
state if the polarization is parallel to the symmetry plane 
and to the other irrep if the polarization is perpendicular 
to the plane. 

To finish the symmetry labels of the initial states form- 
ing FS have to be determined. As was discussed the 
movements of the atoms from the tetr. sites are very 
small in the copper-oxide layers and it can be assumed 
that the shape of the initial state is very close to the 
tetr. case, which can also be shown in the band struc- 
ture calculations. Thus the initial state at FS belongs to 
the same irrep as a Bloch state of copper d x 2_ y 2 orbitals. 
A Bloch sum of a Copper d x 2_ y 2 orbitals at a position 
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TABLE II: Dipole selection rules in the planes of high sym- 
metry. 



k f in the mirror plane E 


Repeated zone of fc/y 


. 1. 


2. 


3. 


4. 


5. ... 


Allowed final State irrep 


. E 


E 


E 


E 


E ... 


Allowed initial state if A _L E . 


. E' 


E' 


E' 


E' 


E' ... 


Allowed initial state if A || E . 


. E 


E 


E 


E 


E ... 




k f in the glide plane A 


Repeated zone of kfu 


. 1. 


2. 


3. 


4. 


5. ... 


Allowed final State irrep 


. A 


A' 


A 


A' 


A ... 


Allowed initial state if A _L A . 


. A' 


A 


A' 


A 


A' ... 


Allowed initial state if A \\ A . 


. A 


A' 


A 


A' 


A ... 





(a) 



(b) 



FIG. 2: Bi2212 FS map. Normal incidence, (a) Polarization 
parallel to x-axis. (b) Polarization parallel to y-axis. 



di = 0.251a + 0.5b - 0.303c is defined as 

4>d x 2_ y 2,CuA r ) = ^e 4k " r -^d x2 _ !;2 (r- di - t v ). (5) 

Tv 

This can be decomposed into components belonging to 
different irreducible representations of the little group of 
k byll 



^£L B 2 ( r ) =\ (<t>d x *_ y z,CuAr) 



e- ik - b / 2 ^,_„ 2 ,c„ 2 (r) 



(6) 



where 0d x2 _ !j2 ,Cu 2 i r ) is a Bloch sum of an identical 
d x 2_ y 2 orbital at d2 = —0.251a + 0b — 0.303c. These 
states can be pictured at T as bonding-antibonding states 
which have the opposite dispersion as a function of 
kj,. </>j A ' is the bonding state with lower energy and 

4>£\ 2 the antibonding state. In the repeated zone 




FIG. 3: Photoelectron spectra along high symmetry lines as 
a function of energy, (a) Along the glide plane, perpendicular 
polarization, (b) Along the glide plane, parallel polarization, 
(c) Along the mirror plane, parallel polarization, (d) Along 
the mirror plane, perpendicular polarization. 



scheme the state 



represents the main band i.e. 



the band of the tctr. structure and the state 



the 



shadow band. 

We have performed first-principles simulations of 
ARUPS spectra in Bi2212 using the one-step model of 
photoemissioniSiii. Fig. |2 shows the calculated FS of 
the orth. Bi2212 in the normal incidence setup. Dou- 
ble headed arrows indicate polarization of the incident 
light. Photon energy was 54 eV. Firstly we would like 
to remark that there seems to be no general correspon- 
dence between the intensitities of the main band and the 
shadow band. The effect of selection rules is seen in high 
symmetry directions S and A. The FS map is probed at 
the fermi level and the initial state belongs to represen- 
tation A' in the glide plane and to the representation E' 
in the mirror plane. In Fig. |2(a)| polarization is parallel 
to the mirror plane and neither the shadow or the main 
is visible in E-direction, wher eas o nly the main band is 
visible in A-direction. In Fig. |2(b) polarization is paral- 
lel to the glide plane and both of the bands are visible 
in E-direction, whereas only the shadow band is visible 
in A-direction. This phenomena has also been shown 
experimentally^. The strange behavior of the shadow 
band follows from the variating irrep of the final state in 
the consequent Brillouin zones. 

Fig. shows photointcnsity as a function of binding 
energy and fen. Calculations were performed within the 
normal incidence setup, photon energy was 40 eV. The 
Fermi function was ignored. Fig. Eta) and Etb) show 
results along the glide plane with polarization perpen- 
dicular to the plane(Fig.(a), main band visible) and par- 
allel to the plane(Fig.(b), shadow band visible). Fig. |3fc) 
and[3Jd) show the corresponding results along the mirror 
plane. Along the glide plane the main and the shadow 
band seem to give continuous intencity as a function of 
fen even when crossing the zone barrier. This can again 
be explained with table [D] When an initial state band 
crosses the zone barrier it can be mapped to the other 
band in the first Brillouin zone and it's irrep, as the irrep 
of the final state, will change. 

Furthermore, there is one surprising consequence of the 
selection rule due to glide plane. As this rule states the fi- 
nal states for the shadow and the main band must belong 
to a different irrep along A-line. Consequently, when pic- 
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turing the band structure, the final state for the shadow 
band cannot be an umklapp of the final state for the main 
band, but on the contrary it must lie in some other avail- 
able band with fixed final state energy. This means that 
the intensities of the bands varie e.g. as a function of 
photon energy in an uncorrelated way, which can also be 
seen as a function of binding energy in Fig|3^a) and^c). 
Actually because of the fact that the width (due to 
and k± dispersion 18 ) of the final state is relatively large, 
both the main and the shadow band can be observed with 
the same single photon energy. 

Available experimental data are consistent with our 
predicted intensities due to structural distortions, but 
we also point out that the intensity variation along the 
glide plane strongly masks any weak sign of antiferro- 
magnetism(AFM) induced spectral features. We have 
modelled the molecular field of the AFM interaction with 
planar ordering, which is known to exist in La2Cu04_jpS, 
by an ad hoc parameter u on copper sites. Our compu- 
tations reveal that, to produce shadow bands in ARUPS 
by AFM, relatively large u has to be used. 

But there is one interesting possiblitity. If both the 
structural and the AFM interactions were present, as 
must be with low doping, the selection rules would change 
again. The correct magnetic group should be determined, 
but because the specification of the magnetic moments 
in Bi2212 is not available and because of limitations of 
our program code, we can only discuss the AFM selec- 
tion rules qualitatively. While maintaining the reflection 
symmetry the ad hoc planar ordering breaks the glide 
symmetry, and, consequently, the derived selection rules 
along the glide plane would be no longer valid. In prac- 
tice, this would mean that, in those high symmetry di- 
rections where the structural distortions forbid intensity, 
there could exist signatures of the FS due to the AFM 
interaction. In Fig. 0]we have calculated the effect of the 
AFM on ARUPS spectra. The magnitude of u was 0.9 
eV, which induces an energy gap of the same magnitude. 
In contrast to Fig. along the glide plane (Fig. 0^ and 
Fig. Et>)) there are no longer strict selection rules for 
the polarization, and the main and the shadow band are, 
though weakly, visible with both polarizations. Along 



the mirror plane [Figs. 0J: and^Ji] the selection rules re- 
main strict. We thus urge experimetalists to focus some 
measurements with high flux and precise polarizations to 
look possible signatures of AFM features. A development 
of the AFM gap with lowered doping and maybe other 
signatures of the AFM may be seen in this way. 

We have shown that structural modifications in Bi2212 
from a body centered tetragonal lattice to a base-centered 
orthorhombic lattice change selection rules in ARUPS 
along the high symmetry lines. Because of the glide 
plane in orthorhombic structure, it turns out that the 
irreducible representation of the final state changes alter- 
natively in the repeated zone scheme of fey space. With 
fixed polarization this yields opposite intensity behaviour 
for the main and the shadow FS in the adjacent Brillouin 
zones. The calculated FS is consistent with available 
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FIG. 4: Photoelectron spectra along the high symmetry lines 
as a function of energy, including AFM-effect, u — 0.9 eV. (a) 
Along the glide plane, perpendicular polarization, (b) Along 
the glide plane, parallel polarization, (c) Along the mirror 
plane, perpendicular polarization, (d) Along the mirror plane, 
parallel polarization. 



experimental data, which is strong evidence against the 
antifcrromagnetic scenario, but a possibility of detecting 
AFM features in low doping is proposed. 
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